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O I Abstract 
!>■ 

The utility of lattice discretization technique is demonstrated for solving nonrelativistic 
quantum scattering problems and specially for the treatment of ultraviolet divergences 
in these problems with some potentials singular at the origin in two and three space 
dimensions. This shows that lattice discretization technique could be a useful tool for the 
numerical solution of scattering problems in general. The approach is illustrated in the 
case of the Dirac delta function potential. 

PACS Numbers 03.80. +r, ll.15.Ha, ll.lO.Gh 

The technique of discretization on lattice (hereafter called lattice technique) |1|, 0, Q has 
been successfully used to deal with ultraviolet divergences in gauge field theoretic problems 
in perturbative expansion, specially those in quantum electrodynamics (QED) and quantum 
chromodynamics (QCD). These ultraviolet divergences in perturbative quantum field theory 
can be eliminated by lattice technique to yield a scale. Except in some simple cases the 
lattice-regularized perturbative series can not be summed up and this makes it difficult to draw 
conclusions about the full solution. The renormalization group (RG) equations 0,0, on the 
other hand, yield many general properties of the full solution from the lattice-regularized results 
of the perturbative expansion. 
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The lattice technique represents a mathematical trick, which removes the ultraviolet diver- 
gences by introducing a cutoff in a regularized Green function. As with any regulator, it is 
removed after renormalization. The physical observables are then obtained in the continuum 
limit, where the lattice spacing is taken to be zero. 

Ultraviolet divergences also appear in the nonrelativistic quantum scattering problems for 
potentials with certain singular behavior at short distances @, ^ § in two and three di- 
mensions. In one dimension these divergences are absent. We show that the application of 
the lattice technique to these potential models leads to a scale and finite physical observables 
after the continuum limit is taken by the usual renormalization procedure. The present work is 
written in a pedagogic style so that it clarifies all the subtleties of lattice technique in a simple 
nonrelativistic problem and should serve as an introduction to the study of lattice technique in 
a complicated field theoretic problem. 

Recently, there have been discussions on renormalization in configuration |^ and momentum 
0, ^, P spaces for potential scattering with the Dirac delta, contact, or zero-range potential. 
In this work the lattice technique is used for potential scattering with delta potential in two 
and three dimensions. In both cases there are ultraviolet divergences. In the two dimensional 
case the divergence is logarithmic in nature, whereas in the three-dimensional case it is linear. 
The lattice-regularized result is finally renormalized and the RG equation written. 

The present potential scattering problem permits analytic solution and is infinitely simpler 
than gauge field theories of QED and QCD where the lattice technique is usually applied. Hence 
the present study will allow us to understand the subtleties of this approach. In gauge field 
theories particles can be spontaneously created and destroyed and the discretization is done 
in the four-dimensional Euclidean space. In potential scattering particle number is conserved 
and we directly discretize the time-independent Schrodinger equation for relative motion in 
three-dimensional Euclidean space. 

The present analytic investigation with delta potential shows the subtleties of the lattice 
technique and demonstrates that this approach can be used for a numerical solution of non- 
relativistic quantum scattering problems in general, not only for two particles but for several 
particles. In this analytic study we calculate the nonrelativistic Green function and the t matrix 
on the lattice. The numerical study remains one to be attempted in the future. 

There is another interest to study the nonrelativistic scattering with delta potential in two 
dimensions. This problem can be considered to be a good model of the ultraviolet structure and 
high energy behavior of field theory 0, ^ 0. Both problems have ultraviolet logarithmic 
divergences, require regularization, are perturbatively renormalizable, collapse for attractive 
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interaction but are asymptotically free, etc. 

We discuss S-wave potential scattering with the delta potential. The partial- wave Lippmann- 
Schwinger equation for the scattering amplitude T{p, q, k'^) in V dimensions at cm. energy k"^ 
is given by 

T{p',p,e) = V{p',p) + J d^qV{p',q)G{q;e)T{q,p,P), (1) 

with the free Green function G{q] k"^) = (/c^ — q"^ + ^0)~^, in units h = 2m = 1, where m is 
the reduced mass. The integral in Eq. (|l]) is over the relevant S'-wave phase space, e.g., we 
take d^q = Airq^dq and d'^q = 27[qdq with q varying from to oo. For the delta potential 
V{p',p) = A, and 

T{p',p,k') = [X-'-Iik)]-\ (2) 

with I{k) = J d^qG{q; k'^). The integral I{k) possesses ultraviolet divergence for V > 1. For 
"D = 3 (2) this divergence is linear (logarithmic) in nature. Finite result for the t matrix of Eq. 
(0) can be obtained only if also diverges in a similar fashion and cancels the divergence of 
lik). 

The solution of the problem can be achieved by discretizing the full Schrodinger equation on 
lattice and finding its solution numerically. Instead, as this problem permits analytic solution, 
we discretize the free Schrodinger equation on lattice and evaluate the lattice-regularized free 
Green function. With the lattice-regularized Green function the ultraviolet divergences are 
avoided. In contrast to the lattice discretization of gauge field theories, where one works in 
terms of Lagrangian densities and path integrals p|, ^, in the present problem it is convenient 
to work in terms of the following time-independent Schrodinger equation for relative motion 

V^0(r) + e4>{r) = 0, (3) 

where the space vector r = {xj),j = The present mathematical treatment is much 

simpler than, but similar to, that in field theory For our purpose we consider the T>- 
dimensional lattice of spacing a. The transition from the continuum to the discrete lattice is 
then effected by making the following substitutions ^ 

Xj — > na = rija, j = 1, 
0(r) ^ 0n = 0(na) 

(j){na + ja) + (j){na — ja) — 2(f){na) . 

Here the space coordinate is discretized by r = na and j is the unit vector in direction j. The 
individual component rij assumes only a finite number of independent values. Outside this 
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range the lattice is assumed to be periodic, so that the nth site can be identified with the 
{n + A^)th site. The active part of the lattice has N'^ sites. 

After discretization, the Schrodinger equation (R|) becomes the matrix equation 



nmV^m 



0, 



(4) 



where 



V 



n+j,m n-j,m 



(5) 



Comparing Eqs. @ and @ we realize that K^m is the discretized version of the operator 
(Vp + k"^)- Hence, the free Green function is the inverse of this operator, defined by 



E-^nm(-^ ^)ml = ^nl- 



(6) 



This inverse operator can be evaluated analytically by working in momentum space where the 
D-dimensional Kroneker 5 functions are represented as 



5. 



V 



(7) 



where q is a dimensionless wave number defined by q = aq with components qj. The integration 
is restricted to the first Brillouin zone — vr < qi < it. In the continuum limit one has the following 
relations for the phase spaces 

1 r 



; .-n/a (27r) 

Using the Fourier representation (^, the matrix K can be written as 



V 



3«q-(n-m) 



V 

E 



-2 



V 

n 



^ dqi 



^iq-(n-m) 



The inverse of the matrix K is now determined by 



[a^k"^ - 2V) + ^2cosgj 



-1) _ „2 
mm 



V 



^ dqi 



□*q-(n-m) 



V 



{a^k^ -2V) + COS qj 



(9) 



(10) 
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This result leads to the following lattice-regularized outgoing-wave Green function 



V 



-2^2 _22?) + ^2cosgj +iO 



:ii) 



With this Green function there is no ultraviolet divergence for a ^ 0. The imaginary part of 
this Green function guarantees unitarity for out going- wave scattering. In the limit a — > 0, the 
regularized Green function reduces to the free Green function: lima^o a; k"^) = G{q; k"^). 

Using the lattice- regularized Green function (pT]), the t matrix (Q) can be rewritten as 



T(/c, A(a), a) 



[\-\a)-lR{k,a)]-\ 



(12) 



where 



lR{k,a) 



(fqGniq, a; k'^ 



V 



'^"-""^ ^ n r 7^%) V^'k^ -2V) + ±2 COS g, + ^0 



1 -1 



(13) 



is a convergent integral for a finite lattice spacing a. In Eq. (|T^) the redundant momentum 
labels p, p' of the t matrix have been suppressed, and the explicit dependences of the t matrix 
on a and A (a) have been introduced. As a — 0, however, this integral develops the original 
ultraviolet divergence. Explicitly, 



\imlR{k,a) 
lim Inlk, a) 



-[c/a + 27r^ik],V = 3, 
27T\n{ak) -in^^V = 2, 



where 



V 

-(n 

1=1 



" dqi 
. (27r)^ 



3 -,-1 

6 — ^2 cos qj , 



(14) 
(15) 

(16) 



is a real finite definite integral. 

Finite results for physical magnitudes, as a — > 0, are obtained from Eq. (|12D if the coupling 
A is also replaced by the so called bare coupling A (a) as in this equation. The bare coupling 
can, for example, be defined by 



\-\a) 



-[c/a + 27r2Ao], P = 3, 
27r[ln(aAo)],P = 2, 



(17) 
(18) 



where Aq is the physical scale of the problem and characterizes the strength of the interaction. 
The quantities A~^(a) of Eqs. (|l^ and (|T^) have the appropriate divergent behavior, as a ^ 0, 
and cancel the divergent part of Iji{k,a) in Eq. ([T^). 
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Regarding the lattice merely as an ultraviolet regulator or cutoff (lattice spacing a) , finally, 
we must take the continuum limit: a — 0. For the present delta potential this limit can be 
taken analytically. In a general problem the limit has to be taken numerically. Both ways 
are illustrated below. After this limit is taken the observables should approach their physical 
values. The question of renormalization is intimately related to the removal of the regulator 
and prediction of physical observables. 

Next the a — > limit is taken analytically in Eq. (|l^) and then we turn to the question of 
renormalization. With the present regularization procedure one has for the lattice-renormalized 
t matrix 

TR{k,\R{A),A) = [\R\A)-XR{k,A)]-\ (19) 

with 

XR{k,A) = \im[lR{k,a)-lRii/A,a)], (20) 

a— »0 

\r\A) = \im[\-\a)-IR{^/A,a)], (21) 

a — >U 

where A is the lattice-renormalization scale of the problem. This scale A should be contrasted 
with the physical scale Aq of Eqs. (|1^ and (pISl). In Eq. (|1^) the explicit dependence of the t 
matrix on both A and the lattice-renormalized coupling Xr{A) has been exhibited. 

After taking the a —>■ limit in Eq. ( pOD and using Eq. (^), the following lattice-renormalized 
function is obtained 



XR{k,A) = lima(2-^)fnr 



dqi 



(27r)3 



-a^A-^ + k^) 



[{a2/c2 -2V + EjLi 2 cos + iO] [-a^A^ -2V + Y.f=i 2 cos q^] 



/A 



A-^ + e 



(A;2-g2 + ,o)(^-2 + ^2)- 



(22) 



Consequently, 



XR{k,A) = -27T^{tk + l/A),V = 3, (23) 

XR{k,A) = 27iln{kA) -in^,V = 2, (24) 

In Eq. (^), if integrals Ir are evaluated and the trivial limit a — > taken, we get 

A^(^) = -[27rV^+27r2Ao]-\P = 3, (25) 

= [27rln(^Ao)]"\l? = 2. (26) 
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The lattice-renormalized coupling for two scales A and are related by the flow equations: 

XR\A) + 27iyA = A^i(A) + 27rVA, (27) 
X^\A) -27r\nA = Xr\Ao) -2n\nAo, (28) 

for V = 3 and 2, respectively. The flow equations are independent of the renormalization 
scheme. 

The present scattering model permits analytic solution and for V = 3 and 2 the exact 
lattice-renormalized t matrices of Eq. (|1^) are given, respectively, by 

Tn{k,\R{A),A) = [\-R\A)+2n\l/A + ik)]-\ (29) 
= [X^\A)-2n\n{kA)+tnY^. (30) 



Explicitly, using definitions ( P^D and (^) for the renormalized coupling, these lattice-renormalized 
t matrices can be written as 

TRik,\niA),A) = [27r\tk-Ao)]-\V = 3 (31) 
= [-27i\n{k/Ao) +in^]-\V = 2. (32) 

These t matrices depend on Xii{A), but not on A, that is the explicit and implicit (through 
\fi{A)) dependences of the t matrix on A cancel. Physics is determined by the value of \r{A) 
at an arbitrary value of ^ |p. For T> =3, the physical scale Aq is related to the scattering 
length ao by = — I/Aq. For V = 2, Aq can also be related to the scattering length [|T^. 

Next we write the RG equation for this problem and show how the limit a — can be taken 
numerically. In this limit the lattice-renormalized t matrix is independent of a, so is invariant 
under the group of transformations a exp(s)a, which form the RG. It is convenient to work 
in terms of the dimensionless coupling, g{a), defined by 

g(a) = c\{a)/a,V = 3, (33) 
= 2n\{a),V = 2. (34) 

The renormalization condition in the a — limit can be expressed as []T| 

a—T(k,g(a),a) = 0, (35) 
da 



or. 



T{k,g{a),a) = 0, (36) 
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where the RG function P{g) is defined by 

dgja) 

m = (37) 

Equation (^) is the RG equation. As the present problem permits analytic solution, the 
constant j3{g) of Eq. (^7|) can be exactly calculated. 

For both V = 3, and 2, (3{g) is a finite quantity independent of a. For V = 3, from Eqs. 
(^), (H) and (0) we have P{g) = -g - g"^. Similarly, for V = 2, from Eqs. (|18D, (H), and 
(13) we have (3{g) = -g^. 

One has the following Taylor series relating the solution for a small finite a, and that for 
a ^ 0: 

T{k, a) = T{k, 0) + ^T"{k, 0) + ^T"'{k, 0)..., 

where prime(s) denote derivative with respect to a at a = 0. Here the linear term in a does not 
contribute, as the RG equation (|35|) yields T'{k, 0) = 0. Though the first order derivative is 
zero by the RG equation, the higher-order derivatives are not zero. Then the converged, a — 0, 
result is given, approximately, by 

T{k, 0) ^ T{k, a) - ^T"{k, a) - ^T"'{k, a)..., (38) 

where the derivatives are to be calculated for a small finite a. For evaluating T" (T'") numeri- 
cally one needs T{k, a) for three (four) adjescent values of a. As more terms are maintained in 
Eq. (|38D a more converged a ^ limit is obtained. 

In summary, we have used the lattice technique for solving the nonrelativistic quantum 
scattering problem with delta potential in two and three dimensions. This technique leads to 
a lattice-regularized Green function. Finite physical result is obtained by employing standard 
renormalization procedure with this regularized Green function as the continuum limit is taken. 
The RG equation is written for this problem. Lattice technique and RG equation should be 
valid for general nonrelativistic potential models with ultraviolet divergence. Though we have 
illustrated the lattice technique for scattering problems with ultraviolet divergences, it should 
be applicable to any scattering problem. In fact, the present study strongly suggests that, 
as in QED and QCD, with the use of modern computers the lattice technique should be a 
powerful alternative tool for the numerical solution of general nonrelativistic few- and many- 
body problems, where, unlike in the present delta potential problems, analytic solutions can 
not be formulated. 

We thank the Conselho Nacional de Desenvolvimento Cientffico e Tecnologico, Fundagao de 
Amparo a Pesquisa do Estado de Sao Paulo, and Financiadora de Estudos e Projetos of Brazil 
for partial financial support. 



8 



References 

[1] Creutz M 1983 Quarks, Gluons, and Lattices, (Cambridge University Press, Cambridge). 

[2] Rothe H 1992 Lattice Gauge Theories: An Introduction, World Scientific Lecture Notes 
in Physics, Vol 43, (World Scientific, Singapore ). 

[3] Creutz M, Jacobs L, and Rebbi C 1983 Phys. Rep. 95C 201 
Rebbi C (1983) Scientific American 248 54 
Kogut J 1983 Rev. Mod. Phys. 55 775 
Callaway D J E 1985 Contemporary Phys. 26 23. 

[4] Wilson K G and Kogut J 1974 Phys. Rep. 12C 78. 

[5] Ryder L H 1985 Quantum Field Theory, (Cambridge University Press, Cambridge). 
[6] Weinberg S (1991) Nucl. Phys. B363 2. 

[7] Adhikari S K and Prederico T 1995 Phys. Rev. Lett. 74 4572. 

[8] Manuel C and Tarrach R 1994 Phys. Lett. B328 113. 

[9] Beg M A B and Furlong R C 1985 Phys. Rev. D 31 1370 
Thorn C 1979 Phys. Rev. D 19 639. 

[10] Adhikari S K 1986 Am. J. Phys 54 362 

Adhikari S K, Gibson W G, and Lim T K 1986 J. Chem. Phys. 85 5580. 



9 



